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We perform density functional theory calculations for the determination of the structural and electronic 
properties of epitaxial graphene on 4H-SiC(0001). Using commensurate supercells that minimize non-physical 
stresses we show that, in contrast with Si-face epitaxial films, the first graphene layer that forms on the C-face 
of SiC is purely metallic with its 7r-bands partially preserved. Typical free-standing characteristics are fully 
recovered with a second graphene layer. We moreover discuss on the magnetic properties of the interface and 
the absence of Fermi-level pinning effects that could allow for a plausible device operation starting from the 
off-state. 



Epitaxial graphene has emerged as a highly attractive 
system for device integration by providing a combina- 
tion of characteristics that constitute a significant ad- 
vantage with respect to competitive graphene produc- 
tion technology: wafer size scaleJi' and direct growth on 
semi-insulating substrates. The process is based on the 
sublimation of Si atoms starting from SiC substrates in 
ultra high vacuum or ambient pressure furnace^. Ac- 
curate control of the growth parameters can allow for 
the reformulation of graphene films from the remaining 
C surface atoms. The structural, electronic and trans- 
port properties of epitaxial graphene strongly depend on 
the the polarization of the SiC surface: Si-face epitaxial 
graphene (i.e. graphene on SiC(OOOl)) is characterized by 
the formation of a first carbon-rich interface layer with a 
6-\/3x 6-y3-R30° surface reconstruction, over which Bernel 
stacked graphene layers grow. Interface interaction im- 
poses significant n-type doping while scatterers reduce 
mobility values with respect to the S'i02-deposited cas^. 
Typical graphene characteristics (e.g. the half-integer 
quantum hall effect) are recovered by the application of 
a gate voltage that lowers the Fermi level around the 
Dirac poinlP. C-face epitaxial graphene (i.e. graphene on 
SiC(OOOl)) is subject to a less stringent rotational order- 
ing with respect to the substratd^ whereas the presence 
of an interface buffer layer, although predicted bythe- 
oretical calculation^^, is still a matter of debat^^'^'li^. 
A complex rotational symmetry of subsequent graphene 
layers other than the AB stacking sequence is presenlPl 
that allows for the manifestation of single-layer proper- 
ties even in the case of a multilayer structur^^i'. C-face 
monolayers show higher mobilities with respect to the 
(0001) case and the typ ical half-integer quantum hall ef- 
fect at low temperature d^ l ^^ l Doping is also present here, 
however controve rsy hol ds over the type and the concen- 
tration of carrierpISIiSl 

In this article we focus on C-face epitaxial graphene by 
means of density functional theory calculations. Starting 
from a basis of lattice commensuration for graphene and 
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FIG. 1. (a) Side view of a (5 x 5) graphene monolayer re- 
laxed on a (4 X 4) 4H-SiC(000l) substrate, (b) Color map 
view of the first epitaxial graphene layer showing the bonding 
characteristics with the substrate, where gradual red to blue 
coloring indicates sp^ to sp^ bonding. Dashed lines show the 
periodically reproduced unit cell. 



the SiC substrate, we show that results can significantly 
differ with respect to non-commensurate models^ ^. Con- 
trary to the Si-fac^l^, we find that the first graphene 
layer preserves a purely metallic character with a com- 
plete absence of a bandgap, while partially maintaining 
its TT-bands. However, interaction with the substrate per- 
turbs the typical Dirac cone that only appears with the 
addition of a second graphene layer. We moreover show 
the presence interface ferromagnetism due to strong ex- 
change interactions and an absence of the Fermi-le vel p in- 
ning effects previously reported for (0001) filmsSHlSl, 

We perform ab initio calculations based on the Density 
Functional Theory (DFT) within the Local Spin Density 
Approximation (LSDA) as implemented in the SIESTA 
computational cod^f^. The 4II-SiC epitaxial graphene 
structures comprise of four bilayers of a (4 x 4) SiC sub- 
strate (passivated with H at the bottom) over which sin- 
gle and double layers of graphene grow, forming a (5 x 5) 
supercell that satisfies lattice commensuration and mini- 
mizes non-physical stresses (Calculated lattice constants: 
asic = 3.076A, agr = 2.46A). Contrary to i?30° re- 
constructions, here there exists no rotational angle for 
the graphene supercell with respect to the SiC substrate. 
Such a configuration is often observed in experiments for 
(0001) samples'*. A basis set of double-C valence (plus 
polarization) orbitals has been used for C (Si) and H, 
along with Troulier-Martins pseudopotentialsi^ for the 
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FIG. 2. Bandstructure of a graphene monolayer on a 4H- 
SiC(OOOl) substrate: (a) spin-restricted configuration, (b) 
majority-spin configuration and (c) minority-spin configura- 
tion. 



modeling of ionic cores. Both basis set sensitivity and 
pseudopotentials have been tested to accurately repro- 
duce the bandstructure of hexagonal SiC polytypes and 
graphene. Sampling of the Brilouin zone takes place by 
a 2 X 2 X 1 Monhorst-Pack grid, whereas bandstructure 
is plotted for 130 nonequivalent A:-points along the closed 
r -J> M isT -J> r path. A mesh cutoff energy of 500 Ry 
has been imposed for real-space integration, while struc- 
tures have been relaxed until forces were less than 0.04 
eV/A. 

At the end of the relaxation process, the first graphene 
layer presents a corrugation with thickness i w 1 A (Fig. 
[ija)) and has structural characteristics that present some 
fundamental differences with respect to the buffer layer of 
the non-commensurate a/3 x \/3R30° reconstruction^ ^. 
A periodic reproduction of the unit cell (Fig. [ijb)) 
shows the presence of small islands covering almost the 
four fifths of the epilayer surface that strongly maintain 
sp^-bonding characteristics and represent the areas of 
the honeycomb lattice that are more distant from the 
substrate {dmax ~ 3 A). Below these islands the sub- 
strate relaxes in an ideal surface reconstruction of the 
C-fac^i^ where also substrate carbon atoms strongly sp'^ 
hybridize. Substrate reconstruction is the origin of the 
energetic stability of this structure: we have calculated 
that the ground state energy proposed by our model is 
significantly lower than the one of the -x/S x \^R30° case. 
The sp^ islands of the epilayer are terminated by carbon 
atoms that strongly bond with the substrate and loose 
the sp^ hybridization for the sp^ one. From a model- 
ing point of view, the increase of the lattice periodicity 
for this supcrcell with respect to the x \/3R30° one 
is in favor of sp^ bonding. Indeed, the biggest part of 
carbon atoms that form the graphene epilayer maintain 
their tt orbitals and only few atoms covalently bond with 
the substrate. It should be moreover noted that in this 
geometrical configuration stresses are minimized, since 
the sp^ bonds of the graphene layer have a mean dis- 
tance of C — Csp2 « 1.43 A whereas interface sp^ ones 
relax at C — CgpS « 1.68 A, which slightly exceeds that of 
pure diamond. Moreover, interesting structural compar- 
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FIG. 3. (a) Total density of states of the majority-spin 
(upper) and the minority-spin (lower) for the graphene/4H- 
SiC(OOOl) system, (b) Projected density of states for the 
majority-spin (upper) and the minority-spin (lower) on the 
carbon atoms of the surface SiC bilayer below the graphene 
layer (lines). The contribution of the 2pz orbitals is also shown 
(lines+symbols) . (c) Local magnetization (p-|- — pi) for the 
previous system, where yellow isosurfaces indicate an excess 
of t electrons. 



isons between the present calculations and the scanning 
tunneling microscopy images of Refs. 5 10 for monolayer 
graphene grown on SiC (0001) (3 x 3) can be made (e.g. 
see figure |ljb) of the present manuscript and figure 2(d) 
of Ref. 10). 

Fig. [2] shows the electronic bandstructure of the first 
graphene layer on 4II-SiC(000i). The key aspect of 
this interaction is that the monolayer preserves a purely 
metallic character with its tt bands prevailing within 
the SiC-substrate bandgap, in consistency with electri- 
cal measurements on C-face monolayerJiSl. This feature 
is in a clear contrast with the carbon-rich layer grown 
on the Si-fac^. However, the interaction with the sub- 
strate, notwithstanding weak, perturbs the Dirac cone 
and the tt bands do not preserve the typical linearity of 
free-standing graphene near the Fermi level. A careful 
view of the spin-polarized bandstructure (Fig. [2jb-c)) 
that corresponds to the most stable electronic configu- 
ration shows some important differences with respect to 
the spin-unpolarized one (Fig. [2]ja)). The presence of 
quasi-flat bands near the Fermi level of this system de- 
rives from SiC surface and is localized on C atoms that 
are not bound to the graphene epilayer. From the spin- 
unpolarized bandstructure one can deduce that these 
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FIG. 4. Bandstructurc of two graphene layers on a 4H- 
SiC(OOOl) substrate for (a) the majority-spin and (b) the 
minority-spin configuration. 

states pin the Fermi level at the interface as similarly 
calculated for Si-face grown graphen^Hl (where this phe- 
nomenon is argued to be one of the main factors for 
the strong n-type doping of that systenl^^'-'^^'). However 
the spin-polarized bandstructurc shows a different pic- 
ture: here the surface states are half-filled and strong ex- 
change interactions make them split above and below the 
Fermi level by ^ 0.7 — O.SeV, leaving the graphene/4H- 
SiC(OOOl) system unpinned (in agreement with previous 
calculations of Ref . [6| . 

The presence of significant exchange interactions 
makes necessary a more careful analysis of magnetism is- 
sues in the interface of these systems. Fig. [3] shows spin- 
up (t) and spin-down (|) density of states (DOS) and 
electronic density p configurations for the graphene/4H- 
SiC(OOOl) system. The total density of states (Fig. [sjja)) 
reveals the presence of a big concentration of spin-up 
states from 0.3eV to leV below the Fermi level, and sim- 
ilarly for spin-down states from OeV to 0.5eV above the 
Fermi level. A careful analysis of the projected density of 
states on the carbon atoms of the first SiC bilayer below 
the graphene layer (Fig. [3f^b)) shows that a major com- 
ponent of these two peaks derives from the those carbon 
atoms that are not covalently bonded with the graphene 
overlayer. Moreover, the contributions of the 2pz orbitals 
are almost exclusive in these peaks (the remaining DOS 
contributed by 2pz orbitals that can be found in the sys- 
tem is localized at the graphene monolayer). The sur- 
face t and I band-splitting gives rise to local magnetiza- 
tion phenomena at the heterostructure's interface (Fig. 
[Sjjc)) where a ferromagnetic order is present. C substrate 
atoms that are not covalently bound with the epilayer 
have magnetic moments that range within 0.6 — 0.67/iB, 
whereas those covalently bound are not magnetic. Mag- 
netism holds only for the first bilayer of the substrate 
and rapidly decreases with distance: the third SiC bilayer 
and the graphene epilayers are not magnetic. Interface 
magnetism however influences the electronic bands of the 
graphene epilayer where a shift of ^ O.leV between ^ and 
I bands can be observed (Fig. ^h) and^c)). 

As in the case of (0001) epitaxy, typical free-standing 
graphene characteristics and the presence of the Dirac 
cone are recovered with the second graphene layer (Fig. 



|4]). This layer interacts weakly with the first epilayer, 
maintaining a mean distance that is similar to that of 
Bernel-stacked graphite and shows a smaller corruga- 
tion with a layer thickness t ^ 0.4 A, in accordance 
with measurement^. It should be again noted here that 
the Fermi level remains unpinned by the interface and is 
purely determined by the position of the Dirac cone. 

To conclude, in this study we have presented ab initio 
electronic structure calculations for monolayer and bi- 
layer epitaxial graphene systems grown on 4II-SiC(000i). 
Contrary to epitaxial graphene on SiC (0001), the first 
graphene layer that grows on the C-face of SiC main- 
tains a purely metallic character and an important pres- 
ence of 77 electrons along with a non-negligible coupling 
with the substrate. The particularity of its geometrical 
configuration consists in a corrugated surface where small 
sp^-bonded islands are terminated by carbon atoms that 
covalently bind to the substrate. Below these islands 
also the surface carbon atoms sp^-hybridize, while their 
TT bands are half-filled and present a ferromagnetic order 
due to the presence of strong electron exchange interac- 
tions. Typical graphene-like characteristics are recovered 
with the addition of a second graphene layer. According 
to this picture C-face epitaxial graphene should main- 
tain some important advantages with respect to the Si- 
face one, like the absence of a Fermi-level pinning effect 
that could allow for a plausible device operation starting 
from an off-state. In this sense, doping effects measured 
in the laboratorj-'^ should likely be attributed to envi- 
ronmental or local disorder factors like the presence of 
impurities, rather than considered as an intrinsic char- 
acteristic of the interface. Methodologically, this study 
argues that approaches that go beyond the x v^^SO" 
reconstructiorP21 are fundamental for the correct model- 
ing of graphene on SiC (0001). 
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